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Abstract 

In this paper, we study discrete-time quantum walks on one-dimensional lattices. We find that 
the coherent dynamics depends on the initial states and coin parameters. For infinite size of 
lattice, wc derive an explicit expression for the return probability, which shows scaling behavior 
P(0,t) t^^ and does not depends on the initial states of the walk. In the long-time limit, the 
probability distribution shows various patterns, depending on the initial states, coin parameters 
and the lattice size. The average mixing time AI^ closes to the limiting probability in linear (size 
of the lattice) for large values of thresholds e. Finally, we introduce another kind of quantum walk 
on infinite or even-numbered size of lattices, and show that the walk is equivalent to the traditional 
quantum walk with symmetrical initial state and coin parameter. 

PACS numbers: 03.67.Lx, 05.40.Fb, 03.65.Db, 03.67.Ca 
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I. INTRODUCTION 



Random walk is related to the diffusion models and is a fundamental topic in discussions 
of Markov processes. Several properties of (classical) random walks, including dispersal 
distributions, first-passage times and encounter rates, have been extensively studied. The 
theory of random walk has been applied to computer science, physics, ecology, economics, 
and a number of other fields as a fundamental model for random processes in time [1-3]. 

Quantum random walk, which is a natural extension of the classical random walk, has at- 
tracted a great deal of attention in the scientific community in recent years. The continuous 
interest in the study of quantum random walk can be partly attributed to its broad applica- 
tions in the field of quantum information and computation [4, 5]. Quantum random walks 
can be used to design highly efficient quantum algorithms for quantum computer [6, 7]. For 
example, Grovers algorithm can be combined with quantum walks in a quantum algorithm 
for glued trees which provides even an exponential speed up over classical methods [7-9]. 
Besides their important applications in quantum computation, quantum walks are also used 
to model the coherent exciton transport in solid state physics [10]. This could be done in the 
framework of the tight-binding model, which is equivalent to the so-called continuous-time 
quantum walk on discrete structures [11, 12]. It is shown that the dramatic non-classical 
behavior of quantum walks can be attributed to quantum coherence, which does not exist 
in the classical random walks. 

In the literature, there are two types of quantum random walks: discrete-time and 
continuous-time quantum walks [13]. The main difference between them is that discrete 
time quantum walk requires a "coin"-which is just any unitary matrix-plus an extra Hilbcrt 
space on which the coin acts, while continuous-time quantum walks do not need this extra 
Hilbert space. Aside from this, these two versions are similar to their classical counterparts. 
Discrete-time quantum walks evolve by the application of a unitary evolution operator at 
discrete time intervals, and continuous-time quantum walks evolve under a (usually time- 
independent) Hamiltonian. Unlike the classical case, the extra Hilbcrt space for discrete-time 
quantum walks means that one cannot obtain the continuous quantum walk from the dis- 
crete walk by taking a limit as the time step goes to zero. This is because discrete time 
quantum walks need an extra Hilbert space, called the "coin" space, and taking the limit 
where the time step goes to zero does not eliminate this Hilbert space. Although there is 
no natural limit to go from the discrete to continuous quantum walks for general graphs, 
Ref. [14] offers a treatment of this limit for the quantum walk on the line, where it is pos- 
sible to meaningfully extract the continuous-time walk as a limit of the discrete-time walk. 
The dynamics of quantum walks of both types has been studied in detail for walks on an 
infinite line-for the continuous-time case in Refs. [11, 15, 16] and for the discrete-time case in 
Refs. [17-20], it has been shown that the properties of discrete and continuous time quantum 
walks are different. 

Here, we focus on discrete-time quantum walks (DTQWs). Previous work have studied 
DTQWs on the line and cycles. The behavior of DTQWs on the line is strikingly different 
from the classical random walks because of quantum interference. The variance o"^ of the 
quantum walk is known to grow quadratically with the number of steps, t, oc t^, compared 
to the linear growth, oc for the classical random walk [13, 19]. Since the cycle (or one- 
dimensional lattice) is a line segment with periodic boundary conditions, the solutions of 
quantum walks on cycles can be simplified greatly on consideration of the Fourier space of 
the particle [21]. For a classical random walk on a ID lattice of size N ^ the mixing time 
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Me converges to the uniform distribution in time 0{—N^lne) [21]. Quantum mechanically, 
the probability oscillates forever and in general do not mix even instantaneously. However, 
by defining a time-averaged probability distribution, the quantum walks can mix to the 
uniform or non-uniform distribution. In Ref. [22], Bednarska et al. have studied the long- 
time limiting probability distribution for a Hadamard walk on ID lattice (cycles). They have 
shown that the Hadamard walk converges to the uniform distribution on odd-numbered size 
of cycles but it converges to a nontrivial distribution on even-numbered ID lattices. Previous 
studies related to DTQWs on ID lattices (cycles) focus on a particular choice of the initial 
state and coin parameter [23-25]. Here, we consider DTQWs on ID lattices for various initial 
states and coin parameters, and discuss the effect of these parameters on the properties of 
the quantum dynamics. 

In this paper, we give a systematic study of DTQWs on ID lattice with various initial 
states and coin parameters. We explore the time evolution of the walk, return probability, 
long-time limiting probability and mixing times, and compare the properties for various 
initial states and coin parameters. The paper is organized as follows: In Sec. II we briefly 
review discrete-time quantum walks on regular graphs. In Sec. HI, we derive analytical 
results for DTQWs on ID lattice and find an exphcit formula for the return probabihty. We 
also do computer simulation to implement DTQWs on ID lattice for various parameters, 
and find that the numerical results accurately agree with our analytical results. In Sec. IV, 
we study the properties of mixing times and discuss the influence of the parameters of the 
walk. In Sec. V, we define another kind of quantum walk on infinite or even-numbered size 
of lattices, and prove that the defined walk is equivalent to the traditional quantum walk 
with symmetrical initial state and coin parameter. Conclusions are given in the last part. 
Sec. VI. 

II. DISCRETE-TIME QUANTUM WALKS 

Discrete-time quantum walk was first introduced by Mayer and Aharonov et al. in 
Ref. [26, 27]. Discrete-time quantum walk takes place in a discrete space of positions, 

with a unitary evolution of coin toss and position shift in discrete time steps. Here, we 
define discrete-time quantum walks (DTQWs) on d-regular graph, which is a regular graph 
each vertex has exactly d edges. 

The discrete-time quantum walk on d-regular graph happens on the coin Hilbert space Tic 
and position Hilbert space Tip, the total Hilbert space is given by ?i = 'Hc®'Hp [13]. If the d- 
regular graph has N vertices we have Tip = : i = 1, 2, A^}, Tic = i = 1,2, d}. 

The coin flip operator C and position shift operator S are applied to the total state in 
H at each time step [13]. The coin flip operation C (only acting on Tic) is the quantum 
equivalent of randomly choosing which way the particle will move, then the position-shift 
operation S moves the particle according to the coin state, transferring the particle into the 
new superposition state in position space. For every vertex, all the outgoing edges is labeled 
as 1,2, ... ,d. The conditional shift operation S moves the particle from u to w if the edge 
{v,w) is labeled by j on v's side [13]: 



The evolution of the system at each step of the walk is governed by the total operator. 




(1) 



(2) 
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where / is the identity operator in Hp. Thus the total state after t steps is given by, 

\m) = u'\m)- (3) 

Finally, we obtain the probability distribution, 

p{x,t) = \{e^,xmm' = E K^i ® (4) 

i=l i=l 

For (i-regular graphs, the coin flip matrix C can be of various forms. The most common 
coins analyzed in the field are the Grover coin and the discrete Fourier transform (DFT) 
coin [13, 28]. It is shown that the choice of coin flip C and initial coin state strongly influence 
the behavior of discrete-time quantum walks. In the next section, we will concentrate 
DTQWs on ID lattice. We will choose a general form of the initial coin state and coin flip 
matrix, and derive analytical results for the walk. 



III. DISCRETE-TIME QUANTUM WALKS ON ONE-DIMENSIONAL LATTICE 

DTQWs on the hue has already been analyzed in detail and the equivalence of all unbiased 
coin operators has been noted by several authors [17-20, 29-31] . Here, we continue to use 
this framework and extend the calculations for DTQWs on ID lattices. 



A. Analytical solutions 

In the following, we restrict our attention to DTQWs on ID lattice. Without loss of 
generality, we consider the one-parameter family of coins. 

The value of p = 1/2 corresponds to the Hadamard coins, which is a balanced coin and 
involves the coin into each direction in Tic with equal probability. Suppose the particle was 
initially {t = 0) localized at vertex xq and the initial coin states distributed in the coin 
subspace, 

\iP{0)) = (V^lei) + ^/l^e"^\e2)) \xo), (6) 

where the two parameters a e [0, 1] and (p e [0, 27r). The position shift operator S has the 
following form [13], 

S={J2\'^- ® + E N + ® Ie2>(e2|. (7) 

i i 

The total states and probability distribution P{x,t) after t steps are determined by 

Eqs. (3) and (4). The solution of the problem can be greatly simplifled in the Fourier space. 
The Fourier transformation of the state in particle Hilbert space can be written as, 

1 ^ 

t)) = -= Y, e^-^'=-/^|V^jv(x, t)),xe {1, 2, TV}. (8) 
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Then the time evolution of the states in the Fourier picture turn into a single difference 
equation, 

\i,r,{kA)) = U{k)\i>N{k,t-l)), (9) 
where U {k) is time evolution operator in the Fourier space, 

The solution of (9) is, 

\^^{k,t)) = iJ{k)'\^^{kM. (11) 

where 0)) is the Fourier transformation of the initial state. To evaluate the pow- 

ers of the propagator f/(/i;), it is convenient to diagonalize U{}i) using its eigenvalues and 
eigenstates. Using En{k) and \qn{k)) to represent the nth eigenvalue and orthonormalized 
eigenvector of the propagator {U{k)\qn{k)) — En{k)\qn{k)) , n — 1,2), Eq. (11) can be written 
as, 

2 

\i^N{k,t)) = J2El{k){q,{k)\i^^ik,0))h{k)). (12) 

i=l 

In the above Equation, Fourier transformation of the initial states is given by |V'iv(^, 0)) = 
;^e^^**^^°/^|Co), where the initial coin state |Co) = \/a|ei) + ^/l — ae*'^|e2). By performing 
the inverse Fourier transformation, we obtain the particle state in position representation 
as follows, 

\^N{x,t)) = ^j:t,e-'''"'''/^\^^{k,t)) 
= ^ELe-^^^'^^-^°^/"'E-=i^K^)(^.WI^o)k,(A;)). ^ ' 

Finally, we get the probability distribution, 

|2 I \/„ +\\\2 



P{x,t) = \{ei\iPN{x,t))Y + \{e2\ipN{x,t))\ 
_ 1 

Ar2 



-ELilELe-'"^'^"""")/^ (14) 



><Y.U^',{mrn\qAk)){<iAk)n)?. 

Substituting eigenvalues Ej{k) and eigenvectors \qj{k)) (See Eq. (Bl) in Appendix B) into 
the above equation, we obtain the probability distribution for DTQWs on ID lattice. We 
have performed numerical implementations which confirm the prediction of Eq. (14). It is 
evident that the probability distribution depends on the initial coin states |Co), eigenvalues 
and eigenstates of U (k). In the following, we will use the above equation to report the time 
evolution of the probability distribution for different initial states and coin parameters. 



B. Time evolution 

In this section, we explore the probability distribution according to Eq. (14). Specifically, 
we consider the following initial coin states (a) |Co) = |ei), (b) |Co) = '^il^i) ^ ks)) and 
(c) |Co) = ^(|ei) ± i\e2)) with p = 3/4, p = 1/2 and p = 1/4. 

Figure 1 shows the probability distribution P{x, t) ait — 20 on ID lattice of size N — 40. 
We note that the initial coin states |Co) give strong influence to evolution of P{x,t). For 
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FIG. 1: (Color online) Probability distribution of DTQWs on ID lattice of size = 40 at t = 20. 
The three rows are for p = 3/4 (row 1), p = 1/2 (row 2) and p = 1/4 (row 3) while the three 
columns correspond to the initial state |Co) = |ei) (column (a)), |Co) = :^(|ei) ± |e2)) (column 
(b)) and \Co) = ;^(ki) ^ ^^2)) (column (c)). The initial node is at xq = N/2 = 20. 

the initial states (a) and (b), P{x,t) is asymmetric at the initial states. On the contrary, for 
initial state (c), P{x,t) is symmetric and displays the same distribution for different initial 
phase (j) = ±7r/2. We also note from the figure that, the coin parameter p does not bias 
the walk; whether P{x,t) symmetric or asymmetric is totally determined by the initial coin 
state I Co). 

Another interesting observation is that the velocities of the two counterpropagating peaks 
is different for different values of p. We find that the peaks spread faster for large values of 
p (Compare the row (l)-(3) in Fig. (1)). This result is consistent with the result in Ref. [32] 
where the positions of the peaks vary linearly with the time steps t. In Ref. [32] , the authors 
show that the peak velocity f is a constant value {v oc ^/p) and differs in sign for the two 
directions. 

It is instructive to consider the extreme values of parameter p. If p = 0, the coin flip 
operation C becomes the Pauli X operation, the two states cross each other going back and 
forth, thereby remaining close to initial excited node. If p = 1, the coin flip operation C 
becomes the Pauli Z operation, the two superposition states |ei) and |e2) move away from 
each other without any diffusion and interference. These two extreme cases are not of much 
importance, but they define the limits of the behavior. Intermediate values of p between 
these extremes show intermediate drifts and quantum interference. 

We also studied the evolution of probability distribution P{x,t) on different lattice size, 
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initial states and coin parameters. The results are analogous to the case we have shown. 
The probability distribution generated by DTQWs consists of two counterpropagating peaks. 
Between the two dominant peaks the probability decays hke while outside the peaks the 
decay is exponential. The probability distribution P{x, t) exhibits symmetric or asymmetric 
characteristics depending on the initial coin states. 



C. Return probabilities 

Now, we consider a particular case of the probability distribution, return probability 
P{x = Xq, t), which means the probability of finding the walker at the initial node. In order 
to study the scaling behavior oi P{x — xo,t), we consider return probability on infinite size 
of lattice. In the continuum limit N — > oo, 6{k) = 27rk/N in the eigenvalues and eigenstates 
(See Eq. (Bl)in Appendix B) become quasicontinuous, the return probability P{x — Xo,t) 
in Eq. (14) can be written as an integral form, 



P{X = Xo,t)\N^oo = ^ Em=l I /-^(Ei=l(emkj(^)) 



x{q,{9)\Co)Ej{9))de\ 



2 



+ /:,(emk2(^))(g2(^)|Co)(-l)V*-(^)(i^P 



(15) 



where Ei — e""^^^\ E2 — —e'^^^'> and uj{9) — arcsin{^/p sin 9) are apphed in the above 
equation. 

In Appendix B, we use the stationary phase approximation (SPA) (See Appendix A) to 
calculate the above integral. We find that the integral is finally simplified as, 

P{x = a;o,t)k^oo = I I ^ (16) 

[0 if t e Odd. 

Equation (16) indicates that the return probability show scaling behavior P{x — 
Xo,t)|Ar^oo ~ We note that the return probabihty does not depend on the parame- 
ters (a and (f) in |Co), see Eq. (6)) of the initial states. Particularly, when p — 1/2, we obtain 
the return probability for Hadamard walks: P{xo,t)\N^oo = 2/(7rt). The parameter p only 
affects the coefficient of the scaling t^^. The scaling behavior of P{xo,t)\N^oo is analogous 
to the return probability of continuous-time quantum walks. For continuous-time quantum 
walks on ID lattice, the return probability is given by 7r{t) — | Jo(2i)p ?» sin^(2i-|-7r/4)/(7ri), 
where J„(x) is the Bessel function of the first kind [33, 34]. Thus, both the return probabil- 
ity of the discrete-time and continuous-time quantum walks show the same scaling behavior 

In order to test the prediction of Eq. (16), Fig. 2 shows the return probability P{x = Xq, t) 
on a ID lattice of size N = 200 with p = 1/4, p = 1/2 and p = 3/4 for the first 100 steps. 
In our calculation, we fix the value of p and try to change the initial states of the walks, and 
find that the return probabilities are exact the same. This confirms our conclusion that the 
return probability is independent on the initial coin states. We also show the theoretical 
predictions of Eq. (16) in Fig. 2, which are in good agreement with the numerical results. 
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FIG. 2: (Color online) Return probability P{x = xo,t) on ID lattice of size N = 200 with p = 1/4 
(triangles), p = 1/2 (dots) and p = 3/4 (squares) in the first 100 steps. The lines show the 
predictions of Eq. (16) for p = 1/4 (dashed line), p = 1/2 (solid line) and p = 3/4 (dotted line), 
respectively. Since P{x = xo,t) equals to at odd- numbered step t, we only plot P{x = xo,t) at 
even- numbered step t in the figure. 



D. Long-time limiting probabilities 

In this section, we consider long time averages of the probability distribution. Generally, 
the time-averaged distribution P{x, T) = ^ Ylt=i P{^i converges to a constant value as 
T — > oo. This value is defined as the long-time limiting probability, 

X{x) = lim P{x,T). (17) 
Substituting the eigenvalues Ei = e~'^^^^^^ and E2 = — e^^^^^^ into Eq. (14), we obtain, 

x(e"-W*(e„|gi(A;))(gi(A;)|Co) 

+ (-l)V-W*(e„|g2(fc))(g2(fc)|Co))p 

_ J_ Y^2 Y^TV -2ni{k-k'){x-xo)/N 

Ar2 Z^m=l Z^fc=l,fc'=l ^ (TO) 

x[e--('=)*(e^|gi(A;))(gi(A;)|Co) ^ ^ 

+(-l)V-W*(e™|g2(fc))(g2(A;)|Co)] 
x[e-(fc')t(Co|gi(fc'))(gi(fc')|e„) 
+(_l)*e— (fc')t(Co|g2(fc'))(g2(A;')|e„^)]. 

Our goal is to calculate the long time averages of the probability distribution. Only two terms 
of the product in Eq. (19) contributes if uj{k) = uj{k'), since limr^oo Z^Li e^''**^'^*^'^)"'^'^'^'))* = 
Suj{k),uj{k') and limy^oo |; ^^^(— l)*e=''*'^'^*^^)+'^*^'^'^)* = 0. Thus, the long-time limiting proba- 
bility xi^) can be simplified as, 

aK-^) — N2 l^m=l l^k=l,k'=l ^uj(k),u{k')^ 

x[(e^|gi(A;))(gi(A;)|Co)(Co|gi(fc'))(gi(A;')|e™) (19) 
H^nr\q2{k)){q2mCo){Co\q2{k')){q^{k')\e^)] 
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FIG. 3: (Color online) Long-time limiting probability x(^) ID lattice of size = 40 for 
p = 3/4 (row 1), p = 1/2 (row 2) and p = 1/4 (row 3) with initial state \Cq) = |ei) (column (a)), 
\Cq) = ± |e2)) (column (b)) and \Cq) = ±i|e2)) (column (c)). The walk starts at 

the initial node xq = N/2 = 20. 



x(x) is dependent on the initial coin states |Co) and coin parameter p (note the eigenstates 
\q{k)) depends on p). Here, we report the long-time limiting probabilities for various initial 
states I Co) and coin parameter p. 

We consider the following initial coin states (a) \Cq) = |ei), (b) |Co) = ± |e2)) and 

(c) |Co) = -^il^i) ±^|e2)) with p = 3/4, p = 1/2 and p = 1/4. Fig. 3 shows the long-time 
limiting probabilities x(^) ID lattice of size = 40. For initial state (c), x(^) is a 
symmetric distribution, i.e., x{^) = x{^') H x — xo = xq — x' . x{x) also displays localization 
on the nodes nearing the initial node xq and the opposite node xq = + iV/2 (mod iV). 
For the same initial state, it seems that small values of p lead to a strong localization than 
large values of p (compare different rows in the same column). For initial state (a), has 
large values at the initial node Xq (or opposite node Xq) and its next node Xq + 1 (or xq + 1). 
These high probabilities are exactly equal, i.e., x{xi^ = xi^o + 1) = xi^o) = xi^o + !)• A 
similar phenomena is also observed for the initial state (b) (See column (b) in Fig 3). x(^) 
has large values at the xq's (or xq's) previous or next nodes, depending on the sign of the 
initial state. Here, ^(^^o) = xi^o + 1) = xi^o) = xi^o + 1) holds only when p = 1/2 (See 
Fig. (3(2b))). The probabilities are equal at nodes x and x + N/2 (mod N) for all the initial 
states and values of p. This behavior is a natural consequence of the periodic boundary 
conditions of the ID lattice. The limiting probability for initial state (c) is symmetric at the 
origin node Xq, i.e., ~ ^o) = xi^o ~ However, this is not true for conditions (b) and 
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(c). This feature is related to the initial coin states. Generally, symmetric initial state leads 
to unbiased limiting probability distribution while asymmetric coin state produces biased 
quantum walks. 

In the above analysis, we have studied the limiting probabilities on lattice of size N — 40. 
For odd values of N, we find that x{^) is a uniform distribution for all the initial states and 
values of p > 0. For even values of A^, we find that the distribution pattern depends on the 
parity of N/2. If N/2 is an even number, x{^) has peaks at the nodes nearing origin node 
and the opposite node. However, if N/2 is an odd number, x{^) has a peak nearing the 
origin node but a minimum nearing the opposite node. The minimum in odd N/2 is virtually 
the mirror of the peak nearing the origin node Xq (See Fig. 1 in Ref. [22]). Nevertheless, the 
situation is different for the extreme case p = 1 and p = 0. If p = 1, two superposition states 
move away from each other without any diffusion and interference, the limiting probability 
is a uniform distribution for all the initial states and values of N. li p — 0, the limiting 
probability x{^) is total determined by the initial states |Co) = {\/a\ei) + y/1 — ae"^|e2)). 
More specifically, the limiting probability x(x) for p = is summarized as. 



It is worth mentioning that the limiting probability distributions between the continuous- 
time and discrete-time quantum walks are different. For continuous-time quantum walks 
(CTQWs), the limiting probability only depends on the parity of N (See Eq.(21) and (22) 
in Ref. [34]). On the contrary, the limiting probability for discrete-time quantum walks 
depends on more ingredients. This is because the coin degrees of freedom of DTQWs offer 
a wider range of controls over the evolution of the walk than the continuous-time quantum 
walk. 

IV. MIXING TIME 

As mentioned in the previous section, the time averaged probability distribution P{x, T) 
converges to the limiting probability x(x) as T — > oo. Here, we study this issue using the 
concept of mixing time. Mixing time represents the rate at which the average probability 
distribution P{x,T) approaches its asymptotic distribution x(a;). The average mixing time 
is defined as follows. 



Fig. 4 (a) shows the time dependence of the variation distance V{T) = \P{x, T)—x{x)\ 
on ID lattices of size = 10, iV = 20 and N = 100. For long times, the variation distance 
V{T) oscillate frequently and decays approximately as 1/T. For odd- numbered size of N, 
the probability mixes to the uniform distribution, we also find a similar behavior of V{T) 
(See V{T) vs T for = 11, = 21 and = 101 in Fig. 4 (b)). 

Fig. 5 shows the dependence of the average mixing time on the lattice size A^ with 
different values of threshold e, for the initial state |Co) — 7f (ki) ± i\^2))- For sufficiently 




(20) 



M,^min{T \ ^ T > t,^\P{x,T) - x{x)\ < e}. 



(21) 
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FIG. 4: (Color online) Variation distance V{T) as a function of time T for ID lattices of size 
= 10 (dotted curve), = 20 (solid curve) and N = 100 (dashed curve). The results are for 
Hadamard quantum walks {p = 1/2) with symmetric initial state |Co) = ;^(|ei) it i|e2)). 

large e, the average mixing time is a linear function of A^. However, for small values of e. 
Me shows wild fluctuation around the linear behavior oc A^. 

We also try to compare the average mixing time for different values of p and initial 
states I Co). We consider quantum walks for the initial states (a), (b), (c) with p = 3/4, 
p = 1/2 and p = 1/4. We find that quantum walk for initial states (c) with p = 3/4 
has a smaller mixing time than the other cases considered here. This may suggest 
that quantum walks with symmetric initial states and large values of p mix to the limiting 
probability distribution fast. We hope this conclusion can be used in constructing efficient 
quantum algorithms. 

V. NEW QUANTUM WALK ON ID LATTICE 

In this section, we introduce another kind of quantum walk on ID lattice. The quantum 
walk is defined on an infinite or even-numbered size of lattice. The walk starts at node xq 
with initial coin state Co|) = ao|ei) + bo\e2), we endow "direction" to the edges in the graph. 
We label each edge a direction (|ei) or |e2)), so that the edges incident on every node can 
be labeled as two different directions and every edge between two nodes has the same label 
at either end. Only ID lattices of even-numbered size satisfy this condition. We illustrate 
this kind of labeling in Fig. 6. The walk is evolved into the superposition of the coin space 
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FIG. 5: (Color online) Dependence of mixing time on the lattice size with e = 0.05, e = 0.1, 
e = 0.2 and e = 0.4. The results are obtained using p = 1/2 and symmetric initial state |Co) = 
7l(|ei)±i|e2)). 

by applying the coin flip operation, then the shift operation S' moves the walker according 
to 

S' = E.eGid^ - ® Iei)(ei| + \x + l){x\ ® |e2)(e2|)] 

+ E.eG.(l^ + ® + \x- l)(x| ® |e2)(e2|)] (22) 



= s, + s. 



2, 



where we use 5*1 and 5*2 to denote the two terms in the above equation. We separate the 
total shift operation S into two different flip operators. Si and 5*2, which are applied to two 
different node group Gi = {x\...,x = Xq — 4:,x = Xq — 2,x = Xo,x = Xq + 2,x = Xq + 4, ...} 
and G2 = {x\...,x = xq — 3,x = Xq — l,x = Xq + l,x = Xq + 3,...}, respectively. We 
implement the above process iteratively to realize a large number of steps of the quantum 
walk. The peculiarity of this walk distinguished from the traditional quantum walk is the 
conditional shift operation S. In the traditional quantum walk, the shift operation S moves 
the walker to the same side of the node, regardless of the position of the walker. However, 
in our defined quantum walk, the shift operation S moves the walker toward different sides, 
depending on the position of walker. 

A natural question is "what's the relationship between the two kinds of quantum walks?" . 
To answer this question, we consider the two quantum walks on the same ID lattice. For 
the sake of simplicity, we compare the wave function of the two quantum walks with dif- 
ferent initial states and coin flip matrixes. Concretely, we consider the traditional quan- 
tum walk with initial coin state \Co) = ao\ei) + bo\e2) and shift matrix C in Eq. (5), and 
our defined quantum walk with initial coin state \Cq) = bo\ei) + ao\e2) and shift matrix 

C' = ("^^ ^ — ] . After t steps, suppose the state of the traditional quantum walk 



at node x is \ip{x,t)) = (A^^tlei) + Bx t\^2))\x), then the state of our defined quantum walk 
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(a) 

4 l^1> l^1> 3 (c) 
l^2> |e2> 

le,) |e2)|ei) |e2) le^) le^} |e,) |e2) le^) |e,) 

— — e e e e e— ■• 




FIG. 6: (Color online) Illustration of the edge labeling on lattice of size = 4 (a) and = 3 
(b). In (a), the edges 1 <-> 2, 3 <-> 4 are labeled as direction |ei), edges 2 <-> 3, 1 ^ 4 are labeled 
as direction |e2). In (b), the edges 1 2, 2 ^ 3 are labeled as direction |ei) and |e2), but edge 
1^3 can't be labeled as a single direction, (c) The labeling of coin direction for the traditional 
quantum walk, where the left side of each node is labeled as |ei) and right side labeled as |e2). (d) 
The labeling of our defined quantum walk. |ei) and |e2) are labeled to each edge and two nodes of 
the edge has the same label. 

\ip'{x,t)) is given by, 

u,'(^ +\\ _ / {Bx,t\(^i) + ^x,t\(^2))\x),x eGi if t G Even, , , 

IV^ [x,t)) - I ^A,^^\e,)-B,^t\e2))\x),x e G2 if t e Odd. ^^"^^ 

The above conclusion can be proved using the method of mathematical induction. For 
t = and t = 1, it is easy to see that the wave functions satisfy the above equation. Now 
suppose for t = Tq (Tq > 1) the conclusion also holds, then we obtain the wave function 
\ip{x,To + 1)) at t = To + 1 according to the iterative relations. 



|^(a;,To + 1)) = (^A^+i,To + VI -P^x+i,To)|ei)|x) 

+ {,/T^A,.^^To - ^/pB.-l,To)\e2)\x) (24) 
= ^x,To+i|ei)|a;) + B^^To+i\e2)\x) 

where we use A^^To+i and Bx,To+i to represent the first two terms. Applying S'{i ® C') to 
\iIj\x,Tq)) and using the iterative relations, we obtain the wave function \iIj\x,Tq + 1)) at 
t = To + 1 

|^/.'(x,To + l)) = 

{(v^^x+i,To + - pBx+i,n)\ei)\x) + 
(-v/r^A^_i,To + v/pS^.-i,To)|e2)|a;),x G G2 if To G Even, (25) 
(y/l - pA^.-i,ro - ^ gx-i,To)|ei)|x) + 
(^Ax+i,To + a/1 - P-Bx+i,To)|ei)|x). X G Gi if To G Odd. 

Comparing Eqs. (24) and (26), we have 

|^'(x,To + l)) = 

^a;,To+i|ei)|x) - 5^,To+i|e2)|a;) x G G2, if e Even, (26) 
Bx,To+i\ei)\x) + A^^To+i\e2)\x), X G Gi, if To G Odd. 
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Therefore, our conclusion is also true for t = Tq + 1. According to the mathematical 
induction, the relation holds for all the time steps. This completes our proof. 

According to the deduction, the two types of quantum walks have the same probability 
distribution. This indicates that our defined quantum walk is equivalent to the traditional 
quantum walk with a symmetrical initial state and coin parameter p. We have performed nu- 
merical implementations to realize the quantum walks defined here, and the results support 
our findings. 



VI. CONCLUSIONS 



In summary, we have studied discrete- time quantum walks on one-dimensional lattices. 
We show that the evolution of the quantum dynamics depends on the initial states and 
coin parameters. For infinite size of lattice, we derive an explicit expression for the return 
probability, which shows scaling behavior P{x = xq, t) ~ and does not depends on the 
initial states of the walk. In the long-time limit, the probability distribution shows various 
patterns, depending on the initial states, coin parameters and the lattice size. The average 
mixing time closes to the limiting probability in linear for large values of thresholds e. 
Finally, we define another kind of quantum walk on infinite or even-numbered size of lattices, 
and find that the walk is equivalent to the traditional quantum walk with symmetric initial 
state and coin parameter. 
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APPENDIX A: THE STATIONARY PHASE APPROXIMATION (SPA) 

Stationary phase approximation (SPA) is an approach for solving integrals analytically 
by evaluating the integrands in regions where they contribute the most [35-37]. This method 
is specifically directed to evaluating oscillatory integrands, where the phase function of the 
integrand is multiplied by a relatively high value. Suppose we want to evaluate the behavior 
of function 7(A) for large A, 

/(A) = ^ 1 gix)e-'f(^^dx. (Al) 

The SPA asserts that the main contribution to this integral comes from those points where 
f{x) is stationary [df{x)/dx = f'{x) = 0]. If there is only one point xq for which f'{xo) = 
and d'^f{x)/dx^\xo = f"{xo) ^ 0, the integral is approximated asymptotically by. 

If there are more than one stationary points satisfy [df{x)/dx = f'{x) = 0], then the integral 
/(A) is approximately given by the sum of the contributions [each being of the form given 
in Eq. (A2)] of all the stationary points [36]. 
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APPENDIX B: CALCULATION OF THE RETURN PROBABILITY USING SPA 



The eigenvalues and eigenstates of U (k) can be written as, 

\ui) = [2 - 2^cos{Lo{k) - e{k))] 

x{VT^P\ei) + hVP + e-('^('=)-'^W)]|e2)}, 
\U2) = [2 + 2^cos{e{k)+uj{k))] 

I \ I \ \U2) 

1^1/ ~ (uiiui)' 1^2/ — (tt2|„2); 

where 6* (A;) = 2'Kk/N and sin a; (A;) = y^sin6'(A;). 

In the continuum hmit N ^ oo, the values of 9k — 2Trk/N are quasicontinuous, then the 
return probability P{x — xo,t) can be written as the integral form in Eq. (15). Now we 
apply SPA to calculate this integral. In Eq. (15), P{x — xo,t)\N-.oo can be written as. 



P{x^xo,t)\N^^ = h + h, (B2) 



where 



+ r(ei|g2(^))(?2(^)|Co)(-l)V*-Wd^r, 
/2 = ^iy (e2|gi(^))(gi(^)|Co)e-^*-Wd^ 
r(e2|g2(^))(g2(^)|Co)(-l)V*-(^)ci^|2. 

J —TT 



+ 

The stationary points of the above integrals satisfy oj'{0) = d axcsinlp sin 6) / d6 = 
■s/pcos{6)/ \/l — psm^{9) = 0. So the contribution of each integral comes from two sta- 
tionary points 9 = 7r/2 and 9 — — 7r/2. The second-order derivations at ^ = ±7r/2 yield 
uj"{±7t/2) = =FVp/V1 - P- According to the SPA and substituting 9 = ±7r/2 into Eqs. (Bl)- 
(B3), we obtain the integral Ii and I2 as follows, 

h ^ \zi + Z2 + Z3 + Z4f, (B4) 

where 

^ ^"-/P ^> (v^ + iVr3^e^'^)e-'*"", 
2^^2^t ^ ' 

= (-1)*-Z2, -^4 = (-l)*'Zl- 

where a;o = arcsin y^. And 

/2~|/l + /2 + /3 + /4r, (B6) 
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If 

is even 



where 

fl = tZl, f2 = -tZ2, h = (-1)V2, U = (-l)Vl. (B7) 

t is odd, zi + Z4 = Z2 + Z3 = and /i + /4 = /2 + /s = 0, the integral equals to 0. If t 
3n, Zi = Z4, Z2 = Zs, fl = f4, f2 = fs, the integral is simplified as, 

I =h + h = |2(zi + ^2)1' + |2(/l + /2)P 

_ V^K . 

Therefore, we obtain the return probability 

p{x = xo,t)k^oo = I if I ^ (B9) 

I if t e Odd. 
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